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We have developed an alternative formulation based on F = Ma rather than F = dP/dt for 
studying variable mass systems. It is shown that F = M a can be particularly useful in this context, 
as illustrated by various examples involving chains and ropes. The method implies the division of 
the whole system into two parts, which are considered separately, allowing to explore certain aspects 
, as constraint forces. 

O 

<N 

on 

£s^j ■ Articles on teaching of variable mass systems have addressed a number of interesting issues 0, 0, 0, 0, 13 ■ F° r 
instance, through adequate examples, Tiersten Q illustrates a sophisticated approach to be applied to open systems, 
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I. INTRODUCTION 



i.e., systems for which there exists a mass influx or efflux. However, in some cases many questions still remain, such 
as the definition of the system and a clarification of the terminology which is used. In addition, the tendency to 
assume that F = dP/dt is, contrary to F = M a, always applicable to a variable mass system is well documented 
and generalized among students. In our opinion these facts are strongly correlated to the way this subject is taught 
^ , in introductory courses. 

& • In this paper we will suggest a different approach, where F = Ma plays a central role. As the whole (closed) 
system under study is separable in two (variable mass) subsystems, we will show that a general equation of the type 
^ ■ F = M a can be used as the equation of motion for one or, in certain cases, for both parts of the system. 
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During the motion, mass is interchanged between both systems although the mass of the total system remains 
constant. Hence, mass is exchanged between the two subsystems and this has to be kept in mind in writing down 
the equations of motion. One may choose the system arbitrarily, including or excluding part of the total number of 
particles, but if interactions act through the system boundary, this must be taken into account. Understanding such 
^) 1 distinctions between different choices of system plays an important role in the subsequent discussion. 

In this context, we develop a method to study variable mass systems which is more clear than the elementary 
traditional method Q . Three sample problems concerning nonrigid systems (ropes and chains) are presented in detail 
\ to illustrate the pedagogical value of the method. 
t-H . We have verified in the classroom that these problems, directly usable in courses of mechanics, provide useful insights 
into the applicability of multi-particle forms of Newton's second law, allowing the teacher to revisit key concepts of 
dynamics. 

The great variety of systems that can be studied by the proposed methodology, such as ropes, chains or conveyor 



, belts, shows that it can be incorporated in the background required in the scientific formation of experts of engineering 



activities. Students should be constantly alert to the various assumptions and approximations in the formulation and 



solution of real problems. Bearing in mind that certain approximations will always be necessary, the ability to 



understand the physical context and to construct the idealized mathematical model for some engineering problems 
is a crucial formative aspect. For instance, the weight of a cable (or a rope) per unit length may be neglected if the 
tension in the cable is much greater than its total weight, whereas the cable weight may not be neglected so far as 
the calculation of the deflection of a suspended cable due to its weight is concerned. The formulation and analysis of 
practical problems involving the principles of dynamics can give the bases to develop systematically these abilities. 

II. GENERAL EQUATIONS OF MOTION 

Newton's second law states that the rate of change of the momentum P of a closed system of mass M is equal to 
the net external force F acting upon it: 
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Alternatively, use can be made of 



M a cm , (2) 



where a cm is the acceleration of the centre of mass of the system. 

If the problem involves variable mass systems, authors assume that it is more convenient to work with momentum, 
and start by using QJ. Of course, under the same conditions, the second statement @ can also be used. 

A. Standard elementary approach 

The traditional method @ to obtain the equation of motion of a variable mass system makes use of Q). The system 
is supposed to have a principal part of mass M moving with velocity v, and a body of mass A M moving with velocity 
u (u > v), which undergo a completely inelastic collision. After the collision the body of mass M + AM is assumed 
to move with a velocity v + A v, so that the variation of linear momentum is (M + AM) (v + Av) — (Mv + AMu). 
Equating the increase in momentum and the impulse of the external force F acting upon the total (constant mass) 
system, and dividing both sides of the equation by At, we obtain, in the limit At —> 0, 

A<Mv, = F + „^. ,3, 

Here, dM/dt is the rate at which the mass enters in the principal part of the system. 

B. Proposed approach and its advantages 

Let us consider a main body which undergoes incremental mass change, increasing (or decreasing) its mass. The 
mass of this body (subsystem II) and its velocity at an arbitrary instant t are M and v, respectively, as it is shown 
in figure 1. The rest of the system (subsystem I) has mass to, moving with velocity u (u > v) at the same instant. 
So, as the motion progresses the two subsystems I and II coalesce to one another. This means that the body I of 
mass m loses mass which immediately becomes part of the main body II. We should keep in mind that to and M are 
time dependent, but the whole system has constant total mass m + M, so that mass is being transferred at the rate 
dM/dt = -dm/dt. 

H V 

M 



FIG. 1: Schematic diagram for the two subsystems I and II. 

Newton's second law for the total closed system reads 

T = (m + M)a cm , (4) 

if we adopt the point of view of 0. T is the external force acting upon the total system. 
The centre of mass velocity is 

' 'mu + Mv) . (5) 



m + M 

The time derivative of this expression (to + M = constant, dm/dt = —dM/dt) is the centre of mass acceleration: 

1 ( du dM dv dAf \ 

m — u + M — + —r- v . (6) 



to + M V'" df df ' df ' df 
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The insertion of this expression in led us directly to 

m , m du AM dv AM , . 

T = F (I) + F (II) = m - — u + M 1 — v , 7 

At At At At w 

where and F^ 11 ) are the forces acting upon subsystems I and II, respectively. 

In order to describe separately the motion of each subsystem I and II, we split J7J) into two: 

F® = TO "^J; for the variable mass system I, (8) 

and 

F*- 11 -* = M— — + (v — u) — — , for the variable mass system II . (9) 
At At 

This splitting of (7J into JSJ and © is physically fruitful as is illustrated in section III. 

The role of the systems I and II can be interchanged, and both of them should be described by the standard 
equation for a variable mass system. With this goal in mind, it is convenient to rewrite JHJ) and (0 using the linear 
momentum. 

From JSJl we obtain 

„m d , . Am .„„, 
F« = - ( tou)-u— , (10) 



and, from JHJ, 



proving the coherence of the treatment. Both expressions have a common generic structure hence, hereafter, the 
symbols I or II are dropped. 

In conclusion, the standard Newton's equation of motion for variable mass systems is obtained: 

dP AM , s 

— = F + u— — . 12 

dt dt v ; 

It is worth pointing out that: 

• the system has instantaneous mass M (to), and linear momentum P — Mv (tou); 

• F is the net external force acting upon the variable mass system; 

• uAM/At (uAm/At) is the rate at which momentum is carried into or away from the system of mass M (to). 

In the present derivation, the definition of the system is clear, as well as the concept of linear momentum which 
agrees with the concept students are acquainted with. Regarding the approach presented previously, we notice that 
F (see ©) is the force acting on the whole (constant mass) system. 

This approach is suitable for undergraduate students and is equivalent to other methods 0, Q • 



III. EXAMPLES 



There are many one-dimensional nonrigid systems, such as chains and ropes, for which the usage of conservation 
laws of energy and linear momentum represent a good approach. However, the conservation laws always refer to a 
definite number of particles, although the system under study might be split into two subsystems of variable mass. 
It is this aspect that is central in our approach. In this context, three illustrative examples are discussed. 
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FIG. 2: The chain with the displacement x of the top end. 



A. Example 1: Falling of a chain 

The upper end of a uniform open-link chain of length I, and mass per unit length X, is released from rest at x — 0. 
The lower end is fixed at point A as it is shown in figure 2. Find the tension T{x) in the chain at point A after the 
upper end of the chain has dropped the distance x. Assume the free fall of the chain. 

Solution. According with the suggested methodology, we divide the chain into two parts each to be treated using 
Newton's second law for variable mass systems i[12|l. For this purpose we consider: 

• subsystem I: the free falling part of the chain, whose mass is X(l — x/2), and velocity v; 

• subsystem II: the rest of the chain, whose mass is Ax/2, and at rest. 
Accordingly, for subsystem I: 



M = X(l - |), P = X(l - |k 



and 



(13) 



System 


M 


P 


F 


u 


uAM/At 


I 


X(l - x/2) 


X(l - x/2)v 


X(l- x/2)g 


V 


-Xv 2 /2 


II 


X x/2 





Xxg/2 + T(x) 


V 


Xv 2 /2 


I+II 


XI 


X(l - x/2)v 


XI g+ T{x) 







TABLE I: Summary of some relevant physical quantities for example 1. See text for notation. 

These quantities are displayed in table I, as well as the corresponding quantities for subsystem II and the total 
system. 

As the net external force is F = X (I — x/2) g, where g stands for the gravity acceleration, equation (|12fl can be 
written as 



(14) 



After straightforward simplifications one concludes that a — dv/dt — g, which is consistent with the free fall 
assumption already used to write F. 

The velocity can then be calculated from dv/dt = v dv/dx = g. The first-order differential equation v dv — gdx = 
is integrated (xq — vq — 0) yielding 



v(x) = {2gxfl 2 . 



(15) 
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As stated above, one can also confirm that a general equation of the type F = Ma applies to subsystem I 
(F = X(l — x/2) g, M = X (I — x/2), a — g). This is a consequence of having v = u. In this case mass is being lost 
but at zero relative velocity. 

To proceed, we apply the reasonings used above to subsystem II. As the net external force on this subsystem is 
F = T(x) + Xxg/2, one finds (see table I), from JTU 



Inserting (|15|1 we obtain for the tension 



= T(x) + l-Xxg + A^-. (16) 



T{x) = ~ i -Xxg. (17) 



One may now check this result using JJJ) for the whole system I+II, and to this end the acceleration of the centre 
of mass must be obtained. 

From figure 2 and table I one has for the centre of mass velocity 

Vom=v(l-^). (18) 

Taking the time derivative we get the centre of mass acceleration: 

= +9(1-%). (19) 



Inserting now v(x) as given by (|15f) in this equation yields 



1.T 

a cm = g{l--). (20) 
Newton's second law (Z), applied to the chain (closed system I+II), reads (see table I) 

q t 

Xlg + T(x) = Xlg(l - — ), (21) 

which confirms (|17(l for T (x). 

This kind of nonrigid systems is also accurate to study energetic processes 0. In fact, the initial mechanical energy 
of the chain is converted into internal energy. 

B. Example 2: Falling of a rope 

A coil of a uniform rope is placed just above a hole in a platform. One end of the rope falls ( without friction or air 
resistance) and pulls down the remaining rope in a steady motion (see figure 3). The rope has total length I and mass 
per unit length X, and starts from rest at x = 0. Find the normal force N(x) exerted upon the coil, and the tension 
T{x) in the rope at a distance x from its lower end. The system is confined in a box so that the rope is prevented from 
lifting up. 

Solution. The rope behaves like a flexible and inextensible one-dimensional conservative system. Under these 
conditions, as it falls down the mechanical energy is conserved. The acceleration of the falling rope is not known, as 
well as the forces N(x) and T(x). Notice that these forces, which result from the constraints imposed on the motion 
of the rope, cannot be specified a priori. 

The acceleration of the falling rope may be obtained form the energy conservation. 

As illustrated in figure 3, the x axis, whose origin is on the platform, points downwards. The variation of the kinetic 
energy of the system is 

AK{x) = ^Xxv 2 , (22) 
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FIG. 3: Scheme of the rope falling off with x being the length of the rope hanging from the hole. 



System 


M 


P 


F 


u 


u AM /At 


I 


X x 


Xxv 


X x g+ T(x) 


V 


Xv 2 


II 


X(l - x) 





X{1 - x)g+ N (x) - T(x) 


V 


-Xv 2 


I+II 


XI 


Xxv 


XI g+ N{x) 







TABLE II: Summary of some relevant physical quantities for example 2. See text for notation. 



where v — dx/dt is the velocity of the moving part of the rope. 

The corresponding variation of the potential energy (whose zero is for x — 0) is 

MJ{x) = -\xgx 2 . (23) 
From the conservation of energy, A E = 0, we obtain the velocity, v, as a function of x, 

v(x) = {gx) 1 ' 2 . (24) 
Taking the time derivative of this equation the acceleration of the falling rope is obtained: 

a dt dx dt 2' 1 ' 

We now illustrate the recommended approach, by looking at the problem from the point of view of a variable mass 
system. 

As figure 3 shows the system at time t can be separated into two subsystems: 

• subsystem I: the hanging part with length x, mass A x and velocity v; 

• subsystem II: the remainder rope at rest on the platform, with mass X(l — x). 

Since we are interested to study the rope as a variable mass system, we put it inside a box, as shown in figure 3. 
Therefore the rope will not depart from the assumed shape as soon as the motion starts. 
For subsystem I (see table II, and figure 3) one has: 

M = Xx, P = Xxv, u — v, and u — — — Xv 2 . (26) 

dt 

Because the external force acting on this system is F — Xxg + T(x), where T{x) is the tension in the rope at the 
boundary of subsystem I, l(T2"|> leads to 

(Xxv) = Xxg + T(x) + Xv 2 , (27) 

and, after using i|25[l. 

T(x) = Xx(a - g) = -^Xxg. (28) 



7 



As we have already said, a general equation of the type F — M a is also valid for this part of the system (F — 
Xxg/2, M = Ax, a = dv/dt — g/2). As in example 1, also here the mass is acquired at zero relative velocity. 

A similar study for subsystem II can be carried out. The relevant physical quantities for this variable mass system 
are also given in table II. 

As the net external force on subsystem II includes the weight of part of the rope, A (7 — x) g, the tension, — T (x) — 
Xxg/2, and the total force by the platform, N(x), the equation of motion Ijl2(l yields 

= X(l - x)g + i Xxg + N (x) - Xv 2 . (29) 

Using (|24|) the normal force on the rope can be obtained: 

N(x)=-X(l-^x)g, (30) 

showing that N(x) = when x — 21/3. We may now confirm this result by applying Newton's second law to the 
closed system I+II using 

The velocity of the centre of mass can be calculated directly from the quantities mentioned in figure 3 and in table 

II: 

"cm = ~v = j {gx) 1 ' 2 . (31) 

The time derivative of this equation, together with (|24|l and <|25|) , allows us calculate the acceleration of the centre 
of mass 

dv cm v 2 x 3 x 



dt = T + l a = 9 2l- (32) 

This equation shows that a cm becomes equal to g when x = 2 1/3, and agrees with N (x) = for that value of x. 
Newton's second law J2J applied to the system I+II (see table II) leads to 

Xlg + N(x) = Xl(^jg), (33) 

resulting a total force on the rope as given by (|30|l . showing the validity of the suggested method. 



C. Example 3: Pulling a rope 

A uniform rope, of length I and mass per unit length X, is pulled along a smooth horizontal surface by a constant 
force F. Find the tension in the rope at any point a distance x from the end where F is applied. 



i " — n 



FIG. 4: Configuration of the rope being pulled on a smooth horizontal surface by the constant force F. 



Solution. Using again the variable mass problem-solving perspective, we present in table III (see also figure 4) 
some of the relevant physical quantities for this example. The origin of the reference system is such that x is the 
coordinate of the point of application of F. 

The equation of motion l|12l) . applied to subsystem I, is written as 

A[A(J - X )v] = T{x) - Xv 2 , (34) 

where T(x) is the tension in the rope at a distance x from the front top of the rope. The terms involving the velocity 
cancel in both sides of (|3"f}l and therefore 



T (x) — X(l — x) a, 



(35) 
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System 


M 


P 


F 


u 


udM/dt 


I 


X(l - x) 


X (I — x) v 


T(x) 


V 


-Xv 2 


II 


X x 


Xxv 


F - T{x) 


V 


X v 2 


I+II 


XI 


Xlv 


F 







TABLE III: Summary of relevant physical quantities for example 3. See text for notation. 



where a is the acceleration of the rope. Regarding subsystem II, wc obtain, from l|12fl 

^ (Xxv) = -T(x) + F + Xv 2 . (36) 

Using l|35[) in this equation the acceleration turns out to be a — F/(Xl). Finally, inserting this acceleration into 
(J2H, we obtain 

T(x)=F(l-^), (37) 

allowing for the correct limits T (x — 0) = F and T (x = I) = 0. 

Since we have now u = v for both subsystems, a general equation of type J5J (see also JSJ) and (0) could be used 
directly leading to the correct results. 

In this context, as already pointed out by Tiersten Q using the concept of flux, equation could not be used to 
obtain the tension T(x). In fact, choosing, for instance, the subsystem I we have: F = T(x), dP/dt = — Xv 2 + 
A (I — x) a, which does not lead to the correct expression for T(x). 



IV. CONCLUSIONS 



In this paper we suggested a method suitable for discussing variable mass systems in a straightforward way and 
presented three applications. Usually the study of the systems under consideration is elaborated over afixed number 
of particles, and so both forms of Newton's second law, either F = Ma or F = dP/dt, can be used 0,0- 

So far as the question of the validity of multi-particle forms of Newton'second law is concerned, three special cases 
must be considered: 

(i) M = constant (in this case F = dP/di and F = Ma are equivalent); 

(ii) u = 0, i.e. the increment of mass d M is at rest; 
and 

(iii) u = v. 

So, as the general equation (|12fl shows, F = dP/diis the correct equation of motion only in two special cases 
(i) and (ii). For example, the raindrop falling through a stationary cloud of droplets : 5] inserts into the case (ii). In 
general, as the present illustrative examples show, conditions (i) and (ii) above are not satisfied and, consequently, 
F = d P/d t does not provide the equation of motion of the corresponding subsystems I and II. 

On the other hand, if the case (iii) applies, a general equation of the type F = Ma can be used (see (0). 
Referring, in particular, to example 3, since in this case v = u for both subsystems, they can be solved by such a 
general equation. The same considerations apply to subsystem I of the other examples. However, let us stress that 
this does not mean that this form of the Newton's second law is more fundamental than in the context of classical 
mechanics, and that care must be taken when referring to variable mass systems. 
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